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O ! 1 INTRODUCTION 



ABSTRACT 

We explore the use of strong lensing by galaxy clusters to constrain the dark energy equa- 
tion of state and its possible time variation. The cores of massive clusters often contain sev- 
eral multiply imaged systems of background galaxies at different redshifts. The locations of 
lensed images can be used to constrain cosmological parameters due to their dependence on 
the ratio of angular diameter distances. We employ Monte-Carlo simulations of cluster lenses, 
including the contribution from substructures, to assess the feasibility of this potentially pow- 
erful technique. At the present, parametric lens models use well motivated scaling relations 
between mass and light to incorporate cluster member galaxies, and do not explicitly model 
line-of-sight structure. Here, we quantify modeling errors due to scatter in the cluster galaxy 
scaling relations and un-modeled line-of-sight halos. These errors are of the order of a few 
arcseconds on average for clusters located at typical redshifts (z ~ 0.2 — 0.3). Using Bayesian 
Markov Chain Monte-Carlo techniques, we show that the inclusion of these modeling errors 
is critical to deriving unbiased constraints on dark energy. However, when the uncertainties 
are properly quantified, we show that constraints competitive with other methods may be ob- 
tained by combining results from a sample of just 10 simulated clusters with 20 families each. 
Cosmography with a set of well studied cluster lenses may provide a powerful complementary 
probe of the dark energy equation of state. Our simulations provide a convenient method of 
quantifying modeling errors and assessing future strong lensing survey strategies. 
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The recent influx of cosmological data from a variety of com- 
plementary observations has led to the development of a highly 
successful concordance cosmology. The overall picture that has 
emerged, where the most significant contributions to the Universe's 
energy content are baryonic matter (~ 5 per cent), non-relativistic 
dark matter (~ 20 per cent), and dark energy (~ 75 per cent), ap- 
pears to be consistent with all current observations. A form of dark 
energy is required within the Friedmann-Robertson-Walker frame- 
work to explain the current epoch of accelerated expansion, which 
is mos t directly evident in the Hubble diagram of Type la super- 
novae dRiess et al.ll998l : |Perlmutter et al.ll999h . However, a variety 
of other combi ned probes such as th e Wilkinson Anisotropy Probe 
(WMAP) (e.g. iKomatsu e"tai]|2010h . baryon acoustic oscillations 
jEfstathiou et alj|2002t ISeliak et al|l2005l : I Eisenstein et al. 1 12001 



baryon fractions ( I Allen et alj|2004h also provide compelling evi- 
dence for its existence. 

Dark energy is typically parameterized by an equation-of-state 
with the form P — w x p, where w x < —1/3 in the current epoch. 
Time variation or observed deviations from w x = — 1 at the present 
day could yield important clues to its nature. In the last decade, 
much effort has been devoted to further constraining the equation of 
state. Owing to the influence of systematic errors in astrophysical 
measurements, the development of complementary observational 
techniques is vital to this task. In this spirit, we investigate the use 
of cluster strong lensing (CSL) as a complementary technique to 
constrain w x . 

The possibility of constraining c osmology with CSL sys- 
tems has been explored in the past (e.g.lPaczvnski & Gorskill 19811: 



cluster abundances jVikhlinin et alj|2009h. cosmic shear measure- 
ments jBacon et alj|2000l: iKaiser et alj|2000l: IVan Waerbeke et~al1 
l2000l : IWittman et al.l I2OO0I : ISemboloni et alj 120061) . and cluster 
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Link & Piercdll998l ;| Cooravlll999l:lGoTse et alj2002l:ISerencl200: 
Soucail et alj |2004 IPaTal et all 120051: IMeneghetti et all 12005 Jbl: 
Macci6ll2005l : iGilmore & Nataraianll2009l : ljullo et alj|2010h . The 
abundance of arcs may provi de usefu l cosm ological constraints. 
IMeneghetti et all d2005bh and iMacciol d2005t) explored the statis- 
tics of arcs in various dark energy cosmologies. They found that 
the variation in arc abundances, particularly at higher redshifts, can 
potentially be used to differentiate between dark energy models. 
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The locations of images in CSL systems also contain useful 
cosmological information. These image positions depend not only 
on the mass distribution, but also on the angular diameter distances 
between the observer, lens, and source. If more than one set of im- 
ages is observed, the geometrical dependence may be exploited to 
probe the cosmological parameters. In this paper we will explore 
the latter technique in greater detail. 

Since the cores of clusters have surface densities which are 
typically much larger than the critical surface density for multiple 
image production, they are prime locations for identifying strongly 
lensed images. Indeed, cluster lenses containing an abundance of 
multiple images have already been utilized to place tight constraints 
on mass distributions in the inner regions of cl usters. The most 
well studied CSL system to date is Abell 1689 l Broadhurst et al.l 
120051 : lHalkola et alj|2006l ; iLimousin et al]|2007l ; ljullo et alJbdToF T 
At present, strongly lensed images from 42 unique sources have 
been identified. Utilizing the large number of constraints and high 
resolution of the Hubble Space Telescope (HST) Advanced Cam- 
era for Surveys (ACS), the parametric mass models cited above can 
reproduce the observed image locations to within a few arcseconds. 



carried out for a growing sample of clusters (e.g. 


Nataraian et al. 


2002 


; Kneibetal. 2003: 


EHasdottir et alj 2007; 


Limousin et al. 


2008 


; Natara 


an et alj 


2009 


; Richard et al. 2009; Smith et al. 2009; 


Richard et al. 


2010dlb 


). 



Preliminary work has been done on bringing these rich sys- 
tems to bear on cosmological parameters. iGolse et al] d2002l) ex- 
plored simulated constraints on the mean matter density O m and 
w x us ing a single cluster with three strongly lensed sources. [Serene! 
j2002l) applied this technique to the cluster C10024+1654 at z = 
0.4. Their res ults were found to be consistent with a flat, acceler- 
ating universe. ISoucail et alj 1 I2OO4T) performed a more detailed ap- 
plication to Abell 2218 using only 4 multiple-image systems. As- 
suming a flat universe, their res ults are consistent with Q, m < 0.3 
wdw x < —0.85. Most recently. I Julio et alj j2oToh obtained Q m = 
0.25±0.05 and w x = — 0.97±0. 07 by combining a strong lensing 
analysis of Abell 1689 with WMAP and X-Ray cluster constraints. 

Theoretical in vestigations of this technique using simulations 
were performed bv lDalal et alj d2005h . They ray-traced through N- 
body simulations to create mock strong lensing image catalogs. 
They then attempted to recover the input cosmology by fitting a 
parametric NFW model and found that significant biases can re- 
sult. Owing to the relatively large errors due to complexities in the 
mass distribution and line-of-sight (LOS) structure, they conclude 
that a single lensing system may not provide significant leverage 
on cosmological parameters. 

Perhaps the most obvious way to overcome these difficulties 
is to obtain a larger sample of CSL systems. This approach has 
the advantage that results obtained from different lines of sight are 
statistically independent. Hence, their results may be combine d in 
a trivial manner. More recently, Gilmore & Nataraiar] j2009h ex- 
plored the prospect of combining CSL systems as a more powerful 
probe of dark energy. They found that competitive constraints can 
be obtained by combining at least 10 lenses with 5 or more image 
systems. Note that space-based measurements, with positional er- 
rors ~ 0.1 arcseconds, and spectroscopic redshift determinations 
(ffz ~ 0.001) are a necessit y. The current work expands upon the 
iGilmore & Nataraianl d2009l) investigation. We further explore the 
issue of overcoming large modeling uncertainties by combining a 
sample of CSL systems. We utilize more realistic simulations of 
cluster lensing systems that include sub-structure in the lens plane 
and line-of-sight structure to address several sources of model- 



ing errors. Moreover, instead of a maximum likelihood routine, a 
Bayesian Markov Chain Monte-Carlo (MCMC) technique is used 
to probe the full parameter space, allowing us to marginalize over 
all lens parameters and draw more robust conclusions on the ability 
of CSL to constrain cosmological parameters. 

This paper is organized as follows. In Section [2] we briefly 
discuss the sensitivity of CSL to cosmological parameters. In Sec- 
tion[3]we describe our simulations and discuss parameter recovery 
through Bayesian MCMC. In Section [4] we use our Monte-Carlo 
simulations to explore modeling uncertainties arising from com- 
plexities in the cluster galaxy population. We also perform a simple 
investigation of errors due to correlated LOS structure, and a more 
detailed analysis of uncorrelated LOS halos. In Section [5] we ex- 
plore the prospect of overcoming these errors using a sample of 10 
simulated clusters. We derive simulated constraints on the dark en- 
ergy equation-of-state for this sample. In Section[6] we explore how 
differences in the large-scale mass distribution of the lens, includ- 
ing bi-modality, affect cosmological constraints. Finally, we offer 
concluding remarks in Section|7] 



2 IMAGE LOCATIONS AND THE DARK ENERGY 
EQUATION-OF-STATE 

In this section, we briefly outline the sensitivity of CS L to cosmo- 
logica l par ameters. For a more detailed discussion, see lGolse et alj 
d2002h and lGilmore & Nataraianl | |2009|) . The lens equation is given 
by 



A = 



D, 



c 2 D ol D . 



■Vcf>{9i 



(1) 



where the angular coordinates of the source i and its correspond- 
ing image(s) are given by j3i and 6i respectively, and cj> is the pro- 
jected Newtonian potential of the lens. For the parametric models 
used in this work, the potential is typically normalized by the as- 
sociated central velocity dispersion, a v . The subscripts o, I, and s 
correspond to the observer, lens, and source. We define D ao to be 
the angular diameter distance from z a to Zb. In the case of a flat, 
two-component universe, D ao is given by 



D(z a , z b ) 



c/Ho 
l + z b 



dz (tt m (l + z) 3 +Q x (z)) 



-1/2 



(2) 



where Ho is the present day Hubble constant, c is the speed of light, 
and f2 m is the present day matter density normalized by the critical 
density. The function Qx (z) is the contribution from dark energy 
and its form depends on the choice of parameterization. In this work 
we consider two parameterizations: 1) a constant equation of state, 
w x . 2) the widely used Chevallier, Polarski, and Linder (CPL) pa- 
rameterization, w x (z) — wo + W a z/(l + z) JChevallier & Polarskil 
1200 ll : lLinderl2003l) . In this case, Q x (z) is given by 



fixO) = n x (i + z 



3(l-\-w +w a ) 



cxp 



3w a z 



1 + 2 



(3) 



Note that the angular diameter distance factor in Q} contains all 
of the explicit dependence on the cosmological parameters, while 
the influence of the mass distribution comes in through the gradi- 
ent term. The CSL technique discussed in this paper is a purely 
geometric probe and utilizes only the angular diameter distances 
for constraints. Aside from the abundance of image systems, the 
influence of cosmology on the detailed properties and formation of 
lensing structures does not come into play. 
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Figure 1. Contours of the ratio rW /r( 2 > in the too — w a plane for lens- 
ing configurations with two multiply-imaged families from background 
galaxies at 2 distinct redshifts. The solid (dashed) curves correspond to 
zi = 0.2(0.3), z 3l = 0.7(1.0), and z s2 = 1.5(4). 



Consider a single source at /3< that is lensed into a set of n 
multiple images 9i,j, where j — 1, n. In this case, it is impos- 
sible to isolate the effect of the cosmological parameters because 
the geometric factor in equation {T} is completely degenerate with 
the normalization of the projected lens potential. If more than one 
source is lensed into multiple image families, as is typically the 
case in CSL systems, then the above degeneracy is broken. Con- 
sider the case of N sources at different redshifts. In this case, there 
is a set of iV lens equations, each given by 



0i = & - Jrr (1) v#fi) 



(4) 



where the T's represent the geometric factors. The superscripts on 
T denote the different source redshifts. Note that equation 10 ap- 
plies to a single cluster lens so that the normalization of <j> is the 
same for each line. Owing to the iV different source redshifts in- 
volved, a change of the cosmological parameters may now be dis- 
cerned from changes in er„. Hence the dark energy equation-of- 
state may be probed using a set of multiply imaged systems as con- 
straints. 

It is instructive to consider the simple example of a lens at zi 
and two families at z B i and z S 2. In this case, the two lens equations 
may be combined by eliminating the normalization of (j>. The effec- 
tive dependence on cosmology for the system is through the ratio 
r (i)/ r (2) Figure □ shows contours of r (1) /T (2) in the w - w a 
plane (CPL model) for t wo different lensing configu rations (for 
the constant w x case, see iGilmore & Nataraiarj|2009h . The solid 
(dashed) curves correspond to z; = 0.2(0.3), z s i = 0.7(1.0) and 
z 3 2 = 1.5(4). The contours, which correspond to one per cent 
differences in the geometric ratios, illustrate degeneracies in the 
wo — w a plane for CSL cosmography with two families. Degen- 
eracies may be broken by using a variety of lensing configurations 
with different lens and source redshifts. 



3 METHODS 

3.1 Simulations of Cluster Strong Lenses 

In this section, we describe how our Monte-Carlo simulations of 
cluster lenses are generated. In the bulk of this work, we em- 
ploy simulated clusters with two components: 1) a smooth large- 
scale density profile. 2) galaxy sized sub-structures within the core 
of the cluster. Lenses with more complicat ed large-scale compo- 
nents are known to exist (for example, see iLimousin et al]|2007l ; 
ljullo & Kneibll2009h . We also address cluster bi-modality in Sec- 
tion!)] 

We model all cluster and galaxy scale halos with a smoothly 
truncated version of the Pseudo-Isothermal Ell iptica l Mass Dis- 
tribut ion (PIEMD) in lKassiola&Kovnerl d 19931) (see lKneib elaTI 
Il99fj) . The profile is characterized by a central velocity dispersion 
a v , core radius r- CO rc, scale radius r cut , and ellipticity parameter 
e. In addition, the profile center (xo,yo) and position angle 6p\ 
are used to characterize lens the orientation. For the circular case, 
r C ore marks the transition from a constant density to p ~ r~ 2 . 
Outside of r cut , the density drops off rapidly as p ~ r~ 4 . The 
truncated PIEMD profile has been used to successf ully model both 
cluster and galaxy scale lenses i n the past (e.g., iNataraian et al.l 
ll998ll2002l ; lLimousin et al.l2007l) . 

The first step in simulating a lensing configuration is to cre- 
ate a smooth, cluster-sized halo. The parameters a v , r core , and 
r cut are drawn uniformly from the intervals 1000 — 1500 km/s, 
30 - 100 kpc, and 800 - 1000 kpc respectively. An ellipticity 
and position angle are randomly assigned in the intervals — 0.3 
and — 360 degrees. The halo center is fixed at the origin. A clus- 
ter redshift in the range z = 0.025 — .6 is drawn from the dis- 
tribution in lGilmore & Nataraianl d2009|). which w as derived from 
the MAssive Cluster Survey dEbeling et alJl200ll) lMACSl. In this 
work, we do not take into account the effects of dark energy on 
cluster assembly characteristics. We note that a change in clus- 
ter proper ties, such as an in c rease in their concentrations as in- 
dicated in iGrossi & SpringeH d2009l) for example, could certainly 
affect the number of multiple-image families and multiplicities ob- 
served. However, we emphasize that the CSL technique utilized 
here is purely geometric in nature, and does not rely on cluster 
properties to constrain dark energy. 

Galaxy-scale potentials are incorporated into the simulated 
clusters with the mean scaling relations derived empirically from 
observations, 



1/4 



Tcut 



1/2 



(5) 



(6) 



where L is the galaxy luminosity. Scaling relations such as these 
are typically used to model galaxy halos in parametric analyses. 
The scaling of the velocity dispersion with luminosity is motivated 
by the Tully-Fisher and Faber- Jackson relations for spiral and el- 
liptical galaxies respectively. We fix the core radii of all galaxy 
sub-halos to be vanishingly small so that the density profiles are 



approximately isothermal (p 



) inside of the scale radii. Al- 



though it is sensitive to the mass of sub-halos, cluster strong lensing 
is at present unable to dis tinguish between different p rofile shapes 
for small-scale structures dNataraian et a l. 2002, 20091) . 

In order to account for cluster-to-cluster variation in the galaxy 
scaling relations, for each lens we draw a* and j-* ut from Gaus- 
sian distributions with mean values a* = 200 km/s, f* ut = 
40 kpc. The standard deviations for the a% and r* ut distributions 
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Figure 2. A simulated lensing configuration with 20 multiply imaged fam- 
ilies. Top panel: a convergence map of the lensing cluster with z\ = 0.24 
and z a = 1.5. The size of the frame is 200 X 200 arcseconds 2 . Bottom 
panel: the image positions, shown as x's, and the critical curves and caus- 
tics for z a = 1.5. The location and redshifts of images from sources lying 
between z s ~ 0.5 — 5 are used to simultaneously invert the lens and con- 
strain the dark energy equation-of-state. 



are 40 km/s and 15 k pc respectively, in acc ordance with observa- 
tional lensing studies jNataraian et alj|2002h . For each galaxy, we 
first draw an J-band absolute magnitude in the range Mi = — 18 
to Mi = -23 from a Schechter function with Mj = -21.7 
and q = —1.3, consistent with recent cluster observations (e.g. 
lHarsono & De ProprisJI 20091) . The PIEMD parameters a v and 
are then calculated using l[5]l and (|6). In section[4] we will draw in- 
dividual a v and r cut values from Gaussian distributions to account 
for scatter in the scaling relations $5$ and {6^. We uniformly draw 
a position angle and an ellipticity in the interval e — — 0.75. 
Finally, a galaxy position is drawn from the cluster-scale PIEMD 



density profile. We place all galaxies within 300 kpc of the clus- 
ter center. In addition to 50 galaxies, for each cluster we draw Mi 
between Mi = —23 and Mi = — 24 from the above Schechter 
function for the Brightest Cluster Galaxy (BCG). For simplicity, 
we use the same scaling relations for the BCG as the other cluster 
galaxies to assign the corresponding a v and r cut . The BCG posi- 
tion is drawn from a uniform distribution within a radius of 50 kpc 
from the cluster center. The BCG a v and r cut parameters are left 
free during parameter recovery. 

The next step is to randomly draw a unique distribution 
of background galaxies to be lens ed by each cluster. We use 
the redshift distribution obtained by iGilmore & Nataraianl d2009h 
from the WFPC2 Hubble Deep Fiel d and Hubble Deep Field- 
South photometric redshift catalogues jFernandez-Soto et al .1 19991: 
lYahata et alj|2000h . We create a population of background galax- 
ies by drawing redshifts and AB(8140) apparent magnitudes in the 
ranges of z = 0.4 - 5.0 and Mab (8140) = 19-28 respectively 
from this distribution. The galaxy positions are drawn uniformly in 
the field of view. The background galaxy number density is fixed 
to 215 arcmin -2 in accord ance with the Hubble Deep Field as in 
IGilmore & Nataraiarj < l2009h . 

After drawing cluster and source properties we lens the pop- 
ulation of background galaxies. We apply a magnitude cut of 
Mab (8140) = 24.5 to remove simulated images that are in prac- 
tice too dim to obtain accurate spectroscopic redshifts. We then 
sort the images into multiple-image families, which consist of at 
least two detectable images. The rest of the images are discarded. 
The positions and redshifts of the final families constitute a fiducial 
catalogue of constraints for the lensing system. In Section [4] we 
will explore perturbations to the catalogs due to complexities in the 
cluster galaxy population and LOS structure. 

Figure [2] shows an example of a lensing configuration gener- 
ated by the above procedure. The top panel shows the convergence 
map of a cluster at zi = 0.24 with z a = 1.5, illustrating the smooth 
cluster component and 51 core cluster galaxies. The field of view is 
200 by 200 arcseconds 2 , consistent with an ACS image. The bot- 
tom panel shows the corresponding critical curves and caustics for 
z s — 1.5. The panel also shows the image positions for 20 strongly 
lensed sources with redshifts ranging from z s — 0.87 to z s = 3.4. 

3.2 Bayesian MCMC 

We use the Bayesian MCMC sampler in the LENSTOOL software 1 
for parameter recovery. In this section we provi de a brief summar y 
of the technique. More details can be found in ljullo et all J2007h . 
LENSTOOL is employed to map image locations to the source 
plane and vice versa using equation (0 for a given lens configu- 
ration. An image plane \ 2 statistic may be computed from 

*' = E E " ^° d ) T c n ~ ^° d ) C7) 

/=! i=l 

where 9 f t and 6^ od are the observed and model positions of image 
i of family / respectively and C is the corresponding covariance 
matrix. Working in the image plane, particularly in the case where 
a large number of \ 2 are computed, is time consuming due to the 
fact that the lens equation must be inverted. A less computation- 
ally expensive approach is to linearize the lens equation and rewrite 
equation (|7) in terms of source plane quantities. In this case, 

1 http://www.oamp.fr/cosmology/lenstool 
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N f n f 

* s = £ £ - (Pf)) T M T CjlM (0 fi - (8) 

/=1 i=l 

where M is the magnification tensor, /3f; is the source plane posi- 
tion, and \Pf \ is the barycenter of family /. In the source plane ap- 
proach, the models with the lowest \ 2 are the those in which source 
plane positions from the same family exhibit minimal scatter. It is 
often the case that some images in a family go undetected due to 
obstruction or demagnification. Note that the location of all images 
in a given family is not required. We only require the identification 
of more than one image in a family. 

The parameters in the models used for the MCMC are sum- 
marized in Table[T] We use flat 50 per cent priors, centered on the 
input values, for all mass profile parameters, except for the clus- 
ter center, which we limit to ±10 arcseconds. Our results are in- 
sensitive to the width of the cluster center limits. For simplicity, 
we assume that mass traces light for the cluster galaxies, and fix 
their positions. In the case of constant equation-of-state, we as- 
sume flat priors of < Q m < 0.7 and —2 < w x < 0. For the 
CPL parameterization, we assume flat priors of —3 < wo < and 
—6 < w a < —wq. We also assume < Q m < 0.7. In practice, 
the exclusio n of w a > — wp woul d come from higher redshift con- 
straints (see lKowalski et al]|2008l) . The CPL parameterization ap- 
proaches wo + w a as z — s> oo. Hence models above the w a — —wo 
line have w x > at high redshifts, implying that dark energy dom- 
inates at early times. Current results favor dark energy making up 
only a small fraction of the total energy density for z S> 1 (e.g. 
iKowalski et alj2008tlKomatsu et allioich . 

The MCMC sampler works by randomly drawing parameter 
sets from Table[T] The \ 2 is computed using the simulated position 
and redshift data and each model is either accepted or rejected us- 
ing a variant of the Metropolis-Hasting algorithm. As the process 
is repeated, the sample of points in the parameter space converges 
to the posterior probability density function (PDF). To avoid get- 
ting stuck in local maxima, the MCMC sampler in LENSTOOL 
util izes the converge nce technique of selective annealing. Follow- 
ing ljullo et alj J200% . we use a convergence rate of 0. 1 and fix the 
number of samples to be 50, 000. 



4 MODELING ERRORS IN CSL 
4.1 Effect of cluster galaxies 

The inclusion of galaxy-scale potentials in cluster lens models 
greatly improves their ability to reproduce observed images. In 
some cases, individual cluster galaxies can have a large impact on 
the configuration of images, and it may difficult to reproduce ob- 
servations without modeling them individually (for example, when 
the presence of a cluster galaxy changes the multiplicity of im- 
ages). However, regarding the bulk of the galaxy population, one is 
forced to make some simplifying assumptions owing to the limited 
number of constraints and large number of parameters required. 
Specifically, the velocity dispersion and scale radii parameters are 
typically assumed to follow empirically motivated scaling relations 
such as (O and ((6). In this section we examine the effect that these 
simplifying assumptions have on our ability to reproduce image lo- 
cations. 

It is well known that the Faber- Jackson and Tully -Fisher rela- 
tions exhibit considerable scatter. We should therefore expect sig- 
nificant scatter in the scaling of a v with luminosity. It is also rea- 
sonable to expect a similar degree of scatter in the r cut relation. The 
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Table 1. Free parameters in the Bayesian MCMC sampler. 



Cluster-scale parameters 


Description 


x,y 


Center of the cluster density profile 


Ov 


Central velocity dispersion 


f core 


Core radius 


rcut 


Scale radius 


e 


Ellipticity 


9PA 


Position angle 


BCG 


&v-> ^cut 


Velocity dispersion and scale radius 




of the brightest cluster galaxy 


Galaxy-scale 




Normalization of the scaling relations 


0} and (|6> 


Cosmological 




Mean matter density 


W x 


Dark Energy equation-of-state 


or 




w ,w a 





deviations of individual cluster galaxies from the relations can in- 
troduce perturbations to nearby images which cannot be accounted 
for with a simple parametric model. We should not expect our mod- 
els to reproduce images to within this perturbation scale. We per- 
form simple Monte-Carlo simulations to quantify these deviations. 

For each of our simulated lensing configurations, we draw 
1000 different galaxy populations. We randomly draw velocity dis- 
persions a v and scale radii r cut from Gaussian distributions with 
mean values obtained through $5$ and ((6) in the fiducial models. 
We use standard deviations of 10, 20, and 30 per cent of the mean 
to explore how varying degrees of scatter affect image positions. 
We lens the given source catalog through each realization and ex- 
amine the scatter in image positions. Figure [3] shows results from 
one such lensing configuration. Each point corresponds to an image 
in one realization of the Monte-Carlo simulations. The open circles 
indicate the position of cluster galaxies. The tangential shapes of 
the image distributions illustrate that the deviations are generally 
larger along the local principal magnification direction. 

In the top panel of Figure [4] we compute image deviations 
relative to the fiducial models with no scatter in the cluster galax- 
ies, and show the probability density of these deviations. The solid, 
dashed and dot-dashed histograms correspond to 10, 20, and 30 
per cent scatter respectively. Since each simulated cluster in our 
sample is at a different redshift, we calculate the deviations in co- 
moving kpc. The average deviation in each case is 2.5, 4.9 and 
6.9 comoving kpc respectively. This corresponds to 0.6, 1.3, and 
1.8 arcseconds for clusters at z ~ 0.2 respectively. In rare cases, 
deviations can be as large as ~ 10 arcseconds. 

We now turn our attention to the effect of the power law in- 
dices in equation ([5) and Although these relations have been 
successfully used in the past for constraining the mass distributions 
of cluster lenses, we examine here more closely the implications 
of these assumptions. As discussed above, the scaling of cr„ with 
~ L 1 / 4 is motivated by well established empirical results. How- 
ever, while a v itself is not directly measured observationally, stud- 
ies of individual galaxy lenses from the SLA CS survey for instance 
(see lBolton et al.ll2008l : lBarnabe et al.l|2009l) suggest that the mea- 
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Figure 3. Results from a Monte-Carlo simulation of CSL illustrating the 
effect of scatter in cluster galaxy scaling relations. We show the lensed im- 
age locations obtained from 1000 random realizations of the cluster galaxy 
population for z\ = 0.34. A 20 per cent scatter in the scaling of velocity 
dispersion and scale radius with luminosity is assumed. The open circles 
show the location of cluster galaxies. The inset (6 X 16 arcseconds 2 ) shows 
the Monte-Carlo realizations for a single image. Each point corresponds to 
the location of the lensed image for a particular realization of the cluster 
galaxy population. 

sured stellar velocity dispersion within an effective radius is a good 
proxy for the velocity dispersion of the lensing mass model. Ow- 
ing to a lack of empirical data, the assumed scaling of r cut with 
luminosity is more uncertain and may vary considerably. 

To determine the effect of assuming an incorrect scaling with 
luminosity, we again perform Monte-Carlo simulations. For each 
cluster in our sample, we draw 1000 realizations where we vary 
the power-law indices. We draw from a Gaussian with mean given 
by 1/4 (a v ) and 1/2 (r cut ). We again use standard deviations of 
10, 20, and 30 per cent of the mean to explore the effect of various 
amounts of scatter. We assume a fixed scatter in the a v and r cu t 
values themselves of 20 per cent. The bottom panel of Figure [4] 
shows image deviations with respect to the fiducial models in which 
there is no scatter in the scalings whatsoever. The 10, 20, and 30 per 
cent cases are shown. We note that varying the power law indices 
has only a mild effect on the deviations if there is scatter present in 
the values themselves. The effect of the indices is essentially lost 
in the scatter. 

Finally, we have checked that the above deviations are rela- 
tively independent of the cosmology, varying only by a few per 
cent for a wide range of Q m and w x parameters. 

4.2 The effect of line-of -sight halos 

4.2.1 The multiple lens plane approximation 

The thin-lens approximation, equation Q}, applies to the case 
where the lensing mass distribution is localized. In this case, the 
deflection angle may be approximated as a single impulse. How- 
ever, in cosmological applications, light from distant sources may 
be deflected by many structures along the path to the observer. The 
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Figure 4. Top panel: magnitude of image deviations, in comoving kpc, 
between models with and without scatter in the cluster galaxy scaling rela- 
tions. For the models with scatter, velocity dispersions and scale radii are 
drawn from Gaussian distributions with mean given by equations j5J and 
{6}. To illustrate the effect of varying degrees of scatter, we assume stan- 
dard deviations of 10, 20, and 30 per cent of the mean values for the solid, 
dashed, and dot-dashed histograms respectively. The scale of these pertur- 
bations is used to quantify modeling errors that arise from complexities in 
the cluster galaxy population. Bottom panel: deviations due to variation in 
the slope of the scaling relations, with a fixed 20 per cent scatter in the 
cr v and r cut values. Assuming mean values of 1 /4 and 1/2 respectively, 
power-law indices for equations {5) and {6} are randomly drawn. The in- 
dices have a minimal effect due to the scatter in a v and r cu t . 



matter distribution may be broken up into multiple lens planes so 
that the deflection angle is approximated by a series of impulses. 
In this case, the angular position of the light ray at lens plane j is 
given by 

i= 1 

where 6\ is the observed image position and a (OA is the deflection 
angle evaluated at the ray position on plane i (see lSchneider et all 
ll992L for example). Note that the lens plane iterator increases with 
distance from the observer. Equation (|9]l may be iterated to trace an 
image location back to the appropriate source plane. 
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Following iHilbert et alj d2009h . we use a quicker and more 
computationally efficient alternative to equation {5J, where the ray 
position on plane j is written in terms of the positions on lens 
planes j — 1 and j — 2, 
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Here, /k is the comoving angular diameter distance, /jj ' = 

D( Za ,Z b )/(l + Z b ). 

We developed our own software to perform the ray tracing 
calculations in the following sections. At present, the code can per- 
form ray tracing through planes consisting of lenses with analytic 
mass profiles, such as the PIEMD and NFW profiles. For the task of 
obtaining image locations, the code traces grid points backwards to 
the appropriate source redshift using equation ( |10t . The set of trans- 
for med grid points in the s ource plane is then divided into triangles 
(c.f. ISchneideret all 19920 . The triangles enclosing the source are 
identified and refined by repeatedly shrinking their sizes to improve 
the resolution of the image location. 



4.2.2 Correlated halos 

In this section, we explore the impact of structures situated along 
the line-of-sight in cluster strong lensing systems. We begin by con- 
sidering, in a simplified manner, the effects of correlated structures. 
In the next section, we will consider uncorrected halos in the light 
cone. 

All CSL studies to date have utilized a single lens plane for 
mass modeling. In what follows, we ask whether correlated struc- 
tures can be reasonably modeled on the same redshift plane as the 
cluster, or whether models consisting of multiple lens planes are 
necessary. We consider the simple example of a cluster lens with 
an infalling group-sized halo residing slightly behind it. Due to 
its proximity, the projection might result in the association of the 
group-sized halo with the cluster system. We simulate a cluster at 
z — 0.192 with an aperture mass of M — 9.5 x 10 14 M s within a 
radius of 1 Mpc. The cluster has 40 core galaxies containing a total 
mass of 2 x 10 13 Mq. The infalling group, with an aperture mass 
of 6.3 x 10 13 Mq, is located 20 comoving Mpc behind the cluster 
at a redshift of z — 0.198, and is situated at (x, y) = (40, 40) arc- 
seconds with respect to the cluster center. We associate 10 galaxies 
with the group. 

We randomly draw 20 different source catalogs with source 
redshifts between z ~ 0.7 — 5 and lens them with the two-plane 
model. Multiply imaged systems are identified and selected. To test 
the effect of using only one lensplane to model the system, we cre- 
ate a model in which the group-sized halo and associated galaxies 
reside on the same lensplane as the cluster. We lens the same 20 
catalogs and compare image locations. 

We find the mean deviation between the two catalogs to be 
0.26 arcseconds, with largest deviations of ~ 1 arcsecond. Such 
cases are rare - 87(97) per cent of deviations are less than 0.5(1) 
arcseconds - and typically occur for images located close to the 
group-sized halo. We also tested a single lensplane model which 
is located halfway between the cluster and group-sized halo. We 




Figure 5. Schematic diagram illustrating the creation of lensplanes to quan- 
tify the effects of LOS halos. A rectangular slice of the Millenium Simula- 
tion box is taken. The locations of halos are projected along the long axis 
and analytic NFW potentials are placed on those sites. The NFW param- 
eters are obtained through scaling relations with mass and redshift. The 
lensplane is inserted at the appropriate redshift and a multi-plane lensing 
algorithm is used to trace rays. In this work, forty-two lensplanes are used 
between 2 = and z = 5. 



found almost no difference in the deviations. Based on the magni- 
tude of deviations, we conclude that in the case where there is a 
dominant lensing system in mass, one can reasonably use a single 
lensplane model. Such an approach may introduce modeling errors 
which are typically on the order of a few tenths of an arcsecond. 
However, we note that these effects are subdominant to the effects 
of scatter in the cluster galaxy population, as shown in the last sec- 
tion. Moreover, in the next section, we will show that the effects of 
uncorrelated galaxy-scale halos along the line-of-sight are greater. 

We now consider the case when the masses of the two line- 
of-sight structures are roughly equivalent. We set the mass of the 
two halos to 4.5 x 10 14 Mq and lens the 20 source catalogs again. 
The deviations increase to 0.91 arcseconds on average and can be 
as high as a few arcseconds. We conclude that in the case with two 
equal-sized halos, it may be necessary to utilize two lensplanes in 
order to minimize modeling errors. 

We note that the above investigation is limited due to the sim- 
plified parametric models used. Galaxy clusters are typically part 
of a rich network of structures which cannot be fully captured in 
the approach we have taken. To test these effects in a more realistic 
way, including larger-scale structures such as filaments, it is best to 
utilize particle data from a cosmological N-body simulation. How- 
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Figure 6. Example convergence map of a simulated line-of-sight created 
from Millenium Simulation halo catalogs, for z s = 1.5. The field of view is 
400 X 400 arcseconds 2 . For each simulation snapshot, a lensplane is created 
by projecting halo locations and placing analytic NFW potentials. The NFW 
parameters are obtained through scaling with the mass and redshift of the 
halo. We use many line-of-sight realizations such as the above to calculate 
perturbations from the cluster-only model, and quantify modeling errors. 



ever, such an investigation is beyond the scope of the current work, 
and we defer to a future paper. 



4.2.3 Uncorrelated halos 

We now turn our attention to the influence of uncorrelated galaxy- 
scale halos along the line-of-sight. Although direct modeling of 
LOS halos from observations may ultimately improve the fit to 
images, perhaps a more practical approach at the present is to 
quantify the errors due to neglecting them in the lens model. We 
utilize halo catalogs from the Millennium Simulation Database 2 
dSprin gel et al. 2005) and the multiple-lensplane approximation 
(see Schneider et al .1 \ 



1992L for example). A schematic of the pro- 



cedure is shown in Figure 5. 

We create one lensplane for each snapshot between z = and 
z = 5, yielding a total of 42 planes. For each lensplane, we take a 
randomly oriented 3-dimensional slice of the corresponding snap- 
shot (see Figure [5] for an illustration). The slices are rectangular 
with length equal to {r CO m(2i-i, z%) + r com (zi, z i+ i)}/2, where 
Zi is the redshift of snapshot i and r com (z a , zt) is the comoving dis- 
tance between z a and Zb. The widths are chosen so that each plane 
fills a 400 x 400 arcsecond 2 area after projection. The locations of 
all halos above a threshold of 10 11 Mq in the slices are projected 
along the long axis. The masses and redshifts in the catalog are used 
to obtain m ean concentration parameters through scaling relations 
obtained bv lGao et alj j200§) . In order to account for scatter about 
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Figure 7. Top panel: ratio of LOS convergences excluding (kfb) and in- 
cluding (k) the contribution from the cluster plane. The contribution from 
foreground and background halos to the total LOS convergence is typically 
less than ~ 10 per cent. Bottom panel: the solid histogram shows image 
deviations, in comoving kpc, between the full light-cone and cluster only 
models. For reference, the dashed histrogram shows the cluster galaxy result 
for 20 per cent scatter in the scaling relations (see Figure |4). Uncorrelated 
LOS halos perturb image locations by ~ 9 kpc on average, corresponding 
to 2.3, 1.6, and 1.2 arcseconds at 2 = 0.2, 0.3 and 0.4 respectively. 



these relations, we draw concentrations fro m a log-normal dis tribu- 
tion with standard deviation given by 0.14 jNeto et alj|2007l) . Cir- 
cularly symmetric NFW potentials with the corresponding param- 
eters are then placed at the location of each halo. Since the mass 
of the NFW profile is not convergent, we truncate each profile at 
the virial radius (for lensing pr operties of the truncated NFW pro- 
file, see iD'Aloisio et alj|2008l) . A typical convergence map for a 
simulated LOS, without cluster, is shown in Figure|6] 

Since the slices are much smaller than the Millennium Sim- 
ulation box, we can perform the above procedure many times for 
each lensplane to create an ensemble. We create random LOS real- 
izations in order to quantify the modeling errors through a Monte- 
Carlo approach. We find that 500 LOS realizations are sufficient 
for numerical convergence. For each lensing configuration in our 
cluster sample, we lens the source catalogs through 500 LOS real- 
izations with the cluster model placed at the appropriate redshift. 
We utilize a ray tracing code that we developed for this purpose. 
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Following lHilbert et al.l J2007b . we calculate the "LOS conver- 
gence" R in order to quantify the lensing effect of the LOS halos. 
Note that this is not equivalent to the multi-plane generalization of 
the convergence, but is nonetheless useful for quantifying the ef- 
fects of LOS structure. We ray trace each image in our simulations 
backwards to their respective source planes. We sum the surface 
mass densities weighted by the appropriate lensing efficiencies at 
each point along the trajectory. We do not subtract off the mean 
density of each plane. For each lensing system, we do this for all 
500 LOS realizations. We calculate the LOS convergence including 
(k) and excluding (kfb) the contribution from the cluster plane. 
The top panel of Figure [7] shows the distribution of kfb /R values. 
The contribution of foreground and background halos to the LOS 
convergence is typically less than ~ 10 per cent. 

The LOS halos perturb the image locations in a manner similar 
to the results shown in Figure [3] though to a larger extent. We note 
that, in some cases, LOS halos can change the multiplicity of im- 
ages. We compare the results of the LOS Monte-Carlo simulations 
to the image catalogs in the cluster-only cases. The solid histogram 
in the bottom panel of Figure [7j corresponds to deviations between 
the full LOS and cluster-only cases. For reference, the dashed his- 
togram shows the case of 20 per cent scatter in the cluster galaxy 
a v and r cut parameters (see the top panel of figure |4j. Figure |7j 
indicates that the uncorrelated halos in the LOS perturb image lo- 
cations by ~ 9 comoving kpc on average (corresponding to 2.3, 
1.6, and 1.2 arcseconds at z = 0.2, 0.3, and 0.4 respectively). 



5 SIMULATED CSL CONSTRAINTS ON DARK ENERGY 

Having explored some of the most significant sources of error in 
CSL, we now use our simulations to perform a feasibility test on 
whether CSL systems with perturbations due to the LOS and clus- 
ter galaxy populations can yield useful dark energy constraints. We 
use the set of ten simulated clusters with 20 families each discussed 
above. Each cluster is placed within an independent light cone gen- 
erated from the Millennium Simulation halo catalogs (see Section 
14.2. 3b - The 50 core galaxies for each cluster are simulated with a 
20 per cent scatter in the luminosity scaling relations (see Section 
I4.lt . Since we have shown that variation in the power-law indices 
of equations l[5} and (|6) makes little difference in the image pertur- 
bations when there is significant scatter in the galaxy population, 
we fix them to the input values. 

We first ray trace the sources through the LOS and cluster 
lensplanes to obtain a catalog of image positions. In order to obtain 
observed positions, we must account for the fact that light from the 
cluster itself is deflected by intervening matter. We therefore lens 
the cluster center through the lensplanes between the cluster and 
the observer, and calculate the image positions with respect to this 
center. We also include the effect of lensing on cluster galaxy posi- 
tions and magnitudes. Finally, image positions are randomly shifted 
by ±2 pixels (0.05 arcseconds per pixel) in order to roughly sim- 
ulate space-based observational errors. We also redraw the source 
redshifts from normal distributions with variance equal to 0.001 to 
simulate measurement errors in spectroscopic redshifts. The final 
catalog of image positions and redshifts of strongly lensed sources 
constitutes the only set of constraints used in parameter recovery. 

We estimate the total errors for each constraint using the 
Monte-Carlo simulations in the last section. Since the Monte-Carlo 
realizations tend to be distributed along the local principal magni- 
fication axis, for each image we transform to the coordinate sys- 
tem in which the magnification tensor is diagonal. In this coor- 
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Figure 8. The effect of underestimating errors in CSL for a single cluster. 
Top panel: the solid contours assume astrometric uncertainties of only 0.1 
arcseconds on the image positions, neglecting all sources of modeling error. 
The dashed contours correspond to parameter recovery using the full error 
estimate, accounting for perturbations due to scatter in the cluster galaxy 
population, LOS halos, and astrometric errors. The input model is depicted 
by a star. Underestimating the total error can lead to erroneous constraints 
on the dark energy equation-of-state. Bottom panel: the solid contours show 
the effect of using photometric redshifts for the high-redshift half of the 
image catalog. For these images, photometric redshift errors are simulated 
by drawing redshifts from a Gaussian distribution with standard deviation 
of 0.5. The full error estimates are used for all images. The dashed contours 
are the same as in the top panel. 



dinate system, the distribution is typically well approximated by 
a bi-variate normal distribution with zero correlations. Note that, 
in addition to the errors, the numerator of each \ 2 contribution is 
calculated in the transformed coordinate system during parameter 
recovery. Modeling errors due to the LOS halos, scatter in cluster 
galaxy scaling relations, and an observational error of 0. 1 arcsec- 
ond are added in quadrature for each image. 
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Figure 9. Simulated CSL constraints on the dark energy equation-of-state 
derived from 10 clusters with 20 multiply-imaged families each. The dark 
contours show the 68, 95, and 99 per cent confidence regions in the ensem- 
ble result. The light contours correspond to the 1 -a regions from individual 
clusters. The stars correspond to the input model (ACDM). The simulated 
image catalogs include perturbations due to scatter in the cluster galaxy 
populations, line-of-sight halos, and observational errors on the image posi- 
tions and redshifts. Top panel: marginalized PDF in the w x and fi m plane. 
Bottom panel: simulated constraints on the CPL parameterization for the 
same ten clusters as the top panel in the w a — wq plane. The line shows 
the border of the w a < — Wq prior. We note that despite the inability to 
reproduce observed image positions to within ~ 1 arcsecond (see Table|2J, 
useful constraints can be obtained by combining results from independent 
clusters. 



Note that we have used the fiducial input models for each clus- 
ter to estimate errors. In practice, one does not have a priori infor- 
mation about the smooth cluster component, cluster-galaxy scaling 
relations, or source positions. This information would be required 
to perform the above Monte Carlo simulations. An initial fit us- 
ing the image positions would therefore be necessary as a starting 
point for the Monte Carlo simulations in this work. Using the ini- 



tial fit, one could then transform the observed images to the source 
plane. Since the initial mass model would not be a perfect fit, im- 
ages from the same source would not necessarily trace back to the 
same source position. One way around this problem is to take the 
barycenters of positions corresponding to single sources. The set 
of barycenters would constitute a model source catalog to base the 
Monte Carlo simulations on. For more details on how Monte Carlo 
based error estimation would be performed in practice, we refer 
the reader to a recent paper on the cluster Abell 1689 that we have 



tne reader to a recent paper on t 
co-authored: ljullo et alj fcOld) . 



We first illustrate the importance of correctly estimating un- 
certainties. Figure [8] shows the contours of the marginalized PDF 
obtained from a single cluster using a mock image catalog of 20 
families which has been perturbed by scatter in the cluster galax- 
ies, LOS halos, and observational errors. From here on, the outer, 
middle and inner contours of PDFs correspond to 99, 95, and 68 
per cent levels respectively. We utilize the simple model described 
in Section[3] whose parameters are listed in Table[T]to perform the 
parameter recovery. The dashed contours show the confidence re- 
gions obtained using the full error estimates, which account for all 
perturbations to image locations in our simulations. The solid con- 
tours in the top panel underestimate the uncertainties by assuming 
only observational errors of 0.1 arcsecond. Note that this is an order 
of magnitude smaller than typical perturbations obtained in Section 
[4] The input model, denoted by a star, is definitively ruled out when 
the errors are underestimated. The bottom panel of Figure [8] shows 
the effect of using photometric redshifts. We simulate errors for the 
high-redshift half of the image catalog by redrawing from Gaussian 
distributions with standard deviations of 0.5. The full positional er- 
rors are used for both sets of contours in the bottom panel. Spectro- 
scopically measured redshifts are necessary to minimize the effect 
of observational uncertainties in CSL cosmography. 

Figure [9] shows the results of combining constraints from our 
simulated ten-cluster sample, where all errors have been fully ac- 
counted for. The top and bottom panels correspond to the con- 
stant equation-of-state and CPL parameterizations respectively. We 
checked the stability of the results by running the MCMC ten times 
for each cluster. Table[2]shows the properties of each cluster in the 
sample. The image plane root-mean-square (RMS) deviations are 
frequently used in strong lensing studies to quantify how well a 
particular model reproduces the observed images. For a given fam- 
ily with n images it is defined to be 



RMS 



(id 



where l o is the observed image position and 6 l is the model image 
position. The total RMS for a given model is obtained by averaging 
i ll lb over all observed families. The last column of Table [2] shows 
that, for a given cluster, the parametric models are not able to re- 
produce image configurations to within ~ 1 arcsecond. In fact, we 
should not expect them to perform any better due to the errors dis- 
cussed in Section|4] However, Figure [9] indicates that the effect of 
these errors can be alleviated by using a relatively small sample 
of clusters. Dark energy constraints from CSL can therefore start 
to be competitive with constraints from other cosmological tech- 
niques upon combining just ~ 10 clusters with ~ 20 families each. 
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Table 2. Properties of simulated clusters. 
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1.19 
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4.8 
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0.70 


0.244 


1.39 
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" The total mass within a circular aperture of 1 Mpc. 
* Averaged over all MCMC samples. Here, we quote the results 
for the constant w x case. 



6 CLUSTER MASS PROFILES 

The main goal of this section is to show that CSL cosmography is 
not limited to a particular type of cluster mass profile. As pointed 
out in Section [2] the CSL technique is a geometric probe that ex- 
ploits the dependence of lensed image locations on angular diam- 
eter distances. In principle, it can be applied to any strong lensing 
system, irrespective of the mass distribution of the lens, as long as 
a sufficient number of constraints are observed to accurately model 
the system. 

Although we have used PIEMD lenses exclusively in this 
work, our results are not contingent upon this choice of profile. 
To illustrate this point, we use the Navarro-Frenk- White profile, 
which has a cusp with the density approaching r _1 in the inner 
regions , and transits to r~ 3 in the outer regions beyond the scale 
radius l lNavarro et al]|l995l Il99rj, Il997l) . This profile is signifi- 
cantly different from the PIEMD case. For lensing properties of the 
NFW profile and a n ell iptical extension of it, w e refer the reader to 
iBartelmand dl996fr > and lGolse & Kneibl J2002t> respectively. 

For comparison, we generate elliptical PIEMD and NFW 
lenses at a redshift of zi = 0.21 with equal masses of 1.6 x 
10 15 Mq. The former has a velocity dispersion of 1300 km/s, core 
radius of 41 kpc and scale radius of 900 kpc. The latter has a con- 
centration parameter c = 4.5 and a scale radius r s — 495 kpc. 
Both have ellipticity parameters of 0.3. Since our purpose in this 
section is solely to compare different cluster-scale mass profiles, 
we neglect the role of sub-structure in what follows. We lens a sim- 
ulated source distribution for both cases and use the same number 
of images as constraints for parameter recovery. We run the MCMC 
sampler using the input models with flat priors. In the NFW case 
all parameters are free. In the PIEMD case, we fix r CO re to the in- 
put value so that the number of free parameters is the same in the 
two cases. Figure [Tol compares the results from both cases. The 



Figure 10. An illustration that similar dark energy constraints may be ob- 
tained from clusters with very different mass profiles. The solid and dashed 
contours correspond to PIEMD and NFW lenses respectively. For both 
cases we use mock catalogs of 21 images. We assume only observational 
errors of 0.1 arcseconds for each image. 
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Figure 11. Same as in Figure [To] except that the dashed contours corre- 
spond to a bimodal cluster consisting of two PIEMD mass peaks. Despite 
having twice the number of free parameters, the bimodal constraints are 
very similar to the unimodal case. 

solid and dashed contours show the PIEMD and NFW constraints 
respectively. We note the similarity of the results, indicating that a 
sample of NFW cluster lenses would generally yield similar con- 
straints to those shown in Section|5] 

In this work we have utilized unimodal lenses consisting of 
only one large-scale halo and a population of galaxy-scale sub- 
halos to model the mass distribution. In reality, cluster lenses can 
be more complex than this, requiring a multi-modal mass distribu- 
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tion to fit the observed images. Here, we address the question of 
whether bimodal clusters (for example, Abell 1689; Abell 2218; 
C10024+16) can be useful for CSL cosmography. 

We use the simulated PIEMD cluster above as a unimodal 
reference. We simulate bimodal lenses with mass equally divided 
between the two mass peaks and total masses equal to the refer- 
ence lens. We generate 10 different lensing configurations for the 
bimodal case using 10 unique source catalogs and random orien- 
tations of the two mass peaks. We then recover the cosmological 
parameters for each realization using the same number of images 
as the reference configuration. Results from the bimodal realiza- 
tion with the widest constraints are shown as the dashed contours 
in Figure QT] The solid contours correspond to the unimodal refer- 
ence lens. Surprisingly, despite the increase in the number of free 
parameters, the bimodal constraints are generally similar to the uni- 
modal constraints. This result is robust for various bimodal config- 
urations and for various ratios of the two mass peaks. The tightness 
of the bimodal constraints is likely due to both the higher amplifica- 
tion of images, on average, and to degeneracies between parameters 
which limit the number of effective free parameters. We conclude 
that despite the larger number of parameters required to model bi- 
modal cluster lenses, these systems may be just as useful for CSL 
cosmography as their unimodal counterparts. Note, however, that 
we have considered a highly idealized model and have not taken 
into account the increased modeling errors that would likely result 
from more complex mass distributions. In addition to the difficulty 
in modeling more complicated systems, the increased errors would 
act to weaken the dark energy constraints and possibly introduce 
biases. We defer a more detailed analysis of more realistic clus- 
ter mass distributions derived from high resolution cosmological 
N-body simulations to future work. 



7 DISCUSSION 

We have shown that, using a relatively small (~ 10 galaxy clus- 
ters) sample of well constrained cluster lenses with parametric mass 
models, it may be possible to obtain constraints on the dark energy 
equation of state and its time variation that are competitive with 
existing probes. Cluster strong lensing may provide a very useful 
complementary tool, particularly in probing dynamic dark energy 
models, which are currently poorly constrained, in the near future. 

We have used simple Monte-Carlo simulations of cluster lens- 
ing configurations to explore some of the potentially largest sources 
of modeling errors. Owing to the large number of parameters that 
would be required to accurately model the cluster galaxy popu- 
lation, parametric models must assume scaling relations between 
galaxy luminosities and their mass profile parameters. While these 
assumptions generally improve agreement with observed images, 
we have shown that scatter in the scaling relations can introduce 
modeling errors as large as ~ 1 arcsecond for clusters at z ~ 
0.2 - 0.3. 

We have used publicly available halo catalogs from the Mil- 
lennium Simulation in order to quantify the lensing effects of inter- 
vening halos along the line-of-sight between the observer and the 
background sources. Owing to limited observational resources at 
the present, most of these halos would not be modeled in practice, 
and could potentially add significantly to the CSL error budget. A 
redshift survey behind the most massive lensing clusters would be 
required to explicitly include structures along the line of sight in the 
modeling. We created an ensemble of line-of-sight realizations con- 
taining analytic galaxy-scale potentials. Using a multi-plane lens- 



ing algorithm, we traced light rays through each line-of-sight real- 
ization, with the simulated cluster placed at the appropriate redshift. 
Although LOS halos typically contribute less than ~ 10 per cent to 
the efficiency-weighted surface mass densities, they can introduce 
deflections with respect to cluster-only models that are typically on 
the order of a few arcseconds for clusters at z ~ 0.2 — 0.3. In rare 
cases, deflections can be as large as ~ 10 arcseconds. 

Given the fact that scatter in the cluster galaxy population and 
LOS halos perturb images typically by a few arc seconds, simple 
parametric models cannot be expected to reproduce image loca- 
tions to within this accuracy. We have shown that underestimating 
the total errors can lead to severe biases in dark energy parameter 
recovery. This is particularly relevant in the case of space-based 
imaging and spectroscopically obtained redshifts, where observa- 
tional errors are typically of the order of ~ 0.1 arcseconds - an 
order of magnitude lower than potential modeling errors. 

We have used our Monte-Carlo simulations to perform a fea- 
sibility test on obtaining competitive dark energy constraints from 
CSL systems. We used a sample of 10 clusters with 20 multi- 
ply imaged sources each. These sources have redshifts between 
z ~ 0.7 — 5. The simulated image catalogs take into account de- 
flections due to scatter in the cluster galaxy scaling relations and 
intervening LOS halos. We also simulated observational errors on 
positions and redshifts. For each simulated cluster lens, we used a 
Bayesian MCMC technique to probe the dark energy equation-of- 
state. We found that the observed images only need to be repro- 
duced to within ~ 1 arcsecond on average to get useful constraints 
upon combining results from ~ 10 clusters. 

Our feasibility test expands on the work of iDalal et all j2005l) 
in a number of ways. Perhaps most importantly, we have shown that 
Bayesian Inference with Markov Chain Monte-Carlo can be used 
to avoid biases that may occur in a best-fit approach by probing 
the full range of parameters that are statistically compatible with 
the image data. Crucial to this point is the use of accurate error 
estimates on image positions, obtained by a detailed analysis of 
each individual cluster, and models containing sufficient degrees of 
freedom to incorporate complexities such as substructure. The large 
range of cosmological parameters compatible with a single lensing 
system can be narrowed by combining results from different sys- 
tems. 

While our simulations have allowed us to explore some of the 
main difficulties in CSL cosmography in a straightforward way, 
there are some important limitations which we plan to address in 
future work. First, our calculation on the effects of the line-of- 
sight includes only mass in collapsed halos within the field of view. 
In reality, the light rays are influenced by an extended network of 
larger scale structures. The mass exterior to halos also contributes 
significantly to the overall deflection of light rays. However, for 
our purposes, the relevant quantities are the relative deflections be- 
tween the observed images, which are typically separated by tens 
of arcseconds. We plan to explore the impact of extended large- 
scale structure in future work exploiting the full particle data and 
periodic boundary conditions of a high resolution cosmological N- 
body simulation. 

Another key limitation is that we have used the same paramet- 
ric profile for the smooth, large-scale component of both the input 
and recovery models. While our focus in this work was primarily 
on the influence of substructure and line-of-sight halos, it will be 
important in future work to test the influence of unc ertainties in 
the lar ge-scale component. In Figure 2 of their paper, IDalal et all 
d2005h show that using a parametric profile different from the input 
model can lead to significant biases in the parameter recovery. We 
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have performed our own tests attempting to recover PIEMD lenses 
with NFW models and vice versa. We found that the Bayesian sam- 
pler did not converge and that the poorness of fit was reflected in 
the evidence values. In practice, it should be straightforward to de- 
termine if the model for the large-scale component generally fails 
to reflect the true mass distribution (i.e. utilizing RMS deviations 
between model and observed images and comparison of Bayesian 
evidence values). However, uncertainties due to less obvious de- 
viations in the large-scale component is a topic that needs to be 
explored. 

In upcoming work, it will also be crucial to test how well para- 
metric models can be used to describe more complex cluster mass 
distributions. In principal, m ore sophisticated m o dels can be intro- 
duced for these systems (e.g. Diego et al .l2005al lblJ Coe et al .120081 ; 
ljullo & Kneibll2009l ; ICoe et alj|20ld) . With ~ 20 families, there 
are certainly enough constraints to increase the number of free pa- 
rameters. However, doing so may degrade the efficacy of derived 
cosmological constraints and introduce additional modeling errors 
and degeneracies. It may be possible to alleviate this problem by 
simply obtaining a larger sample of cluster lenses or including ad- 
ditional independent constraints. Note that we have investigated the 
use of strongly lensed image positions and redshifts, but have not 
considered other possible sources of information. Inclusion of X- 
ray gas temperature profiles, measured cluster velocity dispersions, 
and the stellar velocity dispersion profiles of bright cluster galaxies 
may enable the tightening of derived cosmological constraints with 
more sophisticated mass models. 

We conclude with a brief discussion on observational strat- 
egy. Clearly, the ideal cluster lensing system for cosmography 
would have a relatively simple mass distribution with a large num- 
ber of observed images. However, in practice, there will likely 
be a trade-off between complex lensing systems that produce a 
larger number of strongly lensed images and simple systems with 
fewer constraints. This is due to the possible enhancement of 
stro ng lensing cross sections from substructure an d asymmetries 
(e.g. Bartelmann et al J 19951: Meneghetti et al.ll2007l ). mergers (e.g. 
Torri et alj 2004 Fedeli et all 20061) and line-of- sight matte r (e 



Wambsganss et alj|2005l; [Hennawi et al"]|200?t IffiTbert et alj[200" 
Puchwein & Hilbertll2009l) . The most dramatic and readily observ- 
able cluster strong lenses are likely to be more complicated to 
model. However, in section 6 we showed that a larger number 
of model parameters does not necessarily correspond to weaker 
dark energy constraints. Therefore, it may not be wise to devote 
all resources towards the simplest lensing systems with the fewest 
sources of uncertainty. The complex systems are probably easier to 
find and may provide equal or more leverage on cosmological pa- 
rameters. These benefits may outweigh the costs of more detailed 
analyses. Complex and simple cluster lensing systems are likely to 
play complementary roles in cosmographic applications. 
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